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Abstract 

The FIBRE-pac (FMOS image-based reduction package) is an IRAF-based reduction tool for the fiber 
multiple-object spectrograph (FMOS) of the Subaru telescope. To reduce FMOS images, a number of 
special techniques are necessary because each image contains about 200 separate spectra with airglow 
emission lines variable in spatial and time domains, and with complicated throughput patterns for the 
airglow masks. In spite of these features, almost all of the reduction processes except for a few steps are 
carried out automatically by scripts in text format making it easy to check the commands step by step. 
Wavelength- and flux-calibrated images together with their noise maps are obtained using this reduction 
package. 
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1. Introduction 

FMOS (Kimura et al. 2010) is the near-infrared fiber 
multi-object spectrograph that has been in operation as 
one of the open-use instruments on the Subaru telescope 
since 2009. It can configure 400 fibers of l'/2 aperture in 
a 30' diameter field of view at the primary focus. The 400 
infrared spectra in two groups are taken by two spectro- 
graphs called IRSl and IRS2 (infrared spectrograph 1 and 
2) in either of two modes: a low-resolution mode with a 
spectral resolution of AA ~ 20 A in the 0.9-1.8 /im wave- 
length range, and a high-resolution mode of AA 5 A in 
one of the quarter wavelength ranges of 0.90-1.14, 1.01- 
1.25, 1.11-1.35, 1.40-1.62, 1.49-1.71, or 1.59-1.80 /im. The 
bright OH-airglow emission lines are masked by the mask 
mirror (Iwamuro et al. 2001) installed in these spectro- 
graphs. 

The individual tools for image reduction in this package 
were developed during engineering and guaranteed-time 
observations since 2008, in conjunction with other oper- 
ation and reduction software. Although these tools were 
not developed as part of the public reduction software of 
FMOS, they are now stable, thus they have earned the 
name of "FIBRE-pac" (FMOS Image-Based REduction 
package). The basic concepts underlying the package are 
as follows: 

1. Most of the reduction is processed by IRAF. 

2. Several complicated steps are processed by original 
software written in C using the cfitsio library (Pence 
1999). 

3. Almost all of the reduction processes are automated 
using script files in text format. 



4. Modification of the text files is done by general 
UNIX commands. 

5. The original 2-dimensional information is kept as far 
as possible throughout the reduction processes. 

These concepts enable easy implementation and open pro- 
cessing without the inconvenience of licensed or black box 
parts and ensure traceable operation with visual confirma- 
tion. The 2-dimensional information has advantages not 
only in filtering out unexpected noise using their small 
sizes but also on detection of faint emission-features. In 
this paper, we describe the reduction process of the FMOS 
images based on its Aug. 9-20 11 version, taking care of 
complex conditions in the infrared, using multiple fibers 
and the OH-suppressed spectrograph. 

2. The images 

The FMOS images are acquired by a uniform interval 
non-destructive readout technique called "Up-the-Ramp 
sampling" (hereafter "ramp sampling" for short). A typi- 
cal exposure of 900 seconds consists of 54 images. After an 
exposure is finished, a final frame (treated as a "raw im- 
age" in the reduction process) is created where the count 
in each pixel is calculated by performing a least squares fit 
to the signal count of 54 images. In addition to suppress- 
ing readout noise, the advantages of the ramp sampling 
are 1) saturated pixels can be estimated from the counts 
prior to saturation, and 2) cosmic-ray events can be de- 
tected as an unexpected jump in the counts and removed 
from the final frame (Fixsen et al. 2000). The detection 
threshold of the cosmic-ray events has been currently set 
to lOcr in an empirical way based on the real FMOS im- 
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Fig. 1. Raw image taken in the IRSl low-resolution mode 
with a ramp readout of 54 times using the HAWAII-2 2kx2k 
detector (corresponding to a 900-s exposure). The left half of 
the image consists of spectra in the J band, while the right 
half is in the H band. The typical FWHM of the each fiber 
spectrum is 5 pixels with a pitch of 10 pixels between the 
spectra. 

ages. For IRSl, the fit and the cosmic-ray rejection are 
done during the ramp sampling, so that the final frame 
is ready as the exposure finishes. For IRS2, however, the 
nonlinear bias variation prevents fitting the slope during 
the exposure. Instead, the ramp fitting is executed af- 
ter the simple background subtraction (cf. §3.2) has been 
performed for all of the images taken during the sampling. 
For example, 54 background-subtracted ramp images are 
prepared prior to fitting. Consequently the nonlinear bias 
component is subtracted together with the background 
photons. 

Figure 1 shows the resulting "raw image" after ramp 
sampling 54 times using IRSl. Raw images for IRS2 do 
not exist for the reasons explained above. The raw images 
have the following components: 

1. Thermal background (-300 e/900 s) 

2. Suppressed OH-airglow (—600 e/900 s on average) 

3. Remaining cosmic ray events (—200 e/event) 

4. Read noise (— 10 e rms/54 ramp readout) 

5. Object (-30 e/900 s for a 20 mag(AB) object) 

6. Dark current (<10 e/900 s) 

7. Bias offset (<10 e/exposure) 

8. Cross-talk between quadrants (—0.15% of the count) 

9. Bad pixels having no efficiency, unusual large dark 
current, or too large noise (—0.2% of the pixels) 

10. On-chip amplifier glow at the corners of the four 
quadrants (—2000 e/54 ramp readout) 

11. Unknown external noise in a stripe pattern (quite 
rare and the amplitude is —100 e/900 s at most) 
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12. Point spread function (PSF) of the fiber of 5 pixels 
with a pitch of 10 pixels between spectra (having 
almost no cross talk of the spectra on either side) 

13. Position of an individual fiber along the slit 

14. Optical distortion 

15. Quantum efficiency and its pixel-to-pixel variation 
("ffat" pattern) 

16. System throughput variation with wavelength 

17. Atmospheric transmittance 

18. Intrinsic absorption and emission of the objects 

The thermal background can be removed in the initial sub- 
traction described in subsection 3.2. The 2nd strongest 
component, residual airglow, can be subtracted by inter- 
polation of the background of other fibers after the opti- 
cal distortion is corrected using a dome-fiat image taken 
in the same observation mode. The remaining cosmic- 
ray and other strong noise features having a size smaller 
than the PSF of a fiber are removed in three subsequent 
stages with different threshold counts. Besides the sci- 
ence frames, the following images are necessary for this 
reduction process. 

Detector-flat image : The homogeneous thermal im- 
age of the black board attached to the entrance 
window of the camera dewar. This image gives the 
pixel-to-pixel quantum efficiency ratio of the detec- 
tor. More than 30 sets of the ramp sampling frames 
are averaged to ensure the high signal-to- noise (S/N) 
ratio. This image is included in the reduction pack- 
age together with the other standard images such as 
a bad pixel mask. 

Bad pixel mask : The distribution of bad pixels in the 
detector selected through the reduction process of 
the detector-ffat image. 

Dome-flat image : The dome-flat spectra taken before 
or after observation. This image is used for mea- 
surements of the distortion parameters as well as 
the relative throughput correction of the scientific 
images. 

Th-Ar spectral image : The emission-line spectra for 
wavelength calibration taken just before the expo- 
sures of dome-fiat frames. 

The details of the reduction processes using these images 
are described in the following section. 

3. Reduction Process 

3.1. Preparatory processing 

Before performing a reduction of the science frames, a 
preparatory process is applied to the dome-fiat and Th- 
Ar spectral images to determine the optical distortion and 
the wavelength calibration of the image. 

First, Flat fielding: The dome-fiat image is divided by 
the detector-fiat image to remove differences in the quan- 
tum efficiency between pixels. 

Second, Bad pixel correction: The registered and tem- 
porally prominent bad pixels, which are picked up by sub- 
traction of 3x3 median filtered image, are replaced with 
an interpolated value from the surrounding pixels. 
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Fig. 3. Correction of the spectral distortion of the combined 
ID image of the dome flat. The x-axis of the combined ID 
image (at top) is converted using x' = x -\- {ax + 6), in which 
a and h are coeflicients of a fourth-order polynomial func- 
tion in y. The images of the airglow masks are consequently 
straightened in the corrected image (at bottom). The thick 
lines indicate the shape of the airglow mask at the left and 
the right edges of the detector. 

Third, Correction of the spatial distortion: The y-cods 
of the image is converted using y' = y -\- {aiy + a2)x'^ + 
{biy -\-b2)x. (Here, the origin is at the center of the im- 
age.) The four parameters are chosen to make the PSF 
amphtude in the image projected along the x-axis max- 
imum (cf. figure 2). In this modification process, the 
dome-flat spectrum from each fiber with 9 pixels in width 
is converted into a parallel line to make a "combined ID 
image" (as in figure 3) that includes the pattern of the 
airglow mask. 

Fourth, Correction of the spectral distortion: The x- 
axis of the combined ID image is converted using x' = 
X + {ax + 6), in which the two parameters are determined 
for each line so as to minimize the shift and the magnifi- 
cation difference in the mask pattern. The 200 resulting a 
and b parameters are fitted with fourth-order polynomials 
in y to remove the local matching error in the pattern. 
This conversion process makes the mask pattern straight 
in the column direction (cf. figure 3), which is necessary 
for better subtraction of the residual airglow lines in the 
science frames. 

These parameters obtained from the dome-flat image 
are also applied to the Th-Ar spectral image to make a 
combined ID image, as shown in figure 4. Although the 
slit image of each emission line and of the airglow masks 
should be parallel, alignment errors cause small differences 
between them. 

Finally, Wavelength calibration: The correlation be- 
tween the observed wavelength and the pixels of each spec- 
trum are determined, as represented by A = px^ + qx'^ + 
rx-\-s. The resulting p, r, and s coefficients are used to 
calculate the corresponding wavelengths in each spectrum 
without changing them, because the x positions of the in- 
dividual fibers in the slit may exhibit some scatter due 
to imperfect fiber alignment along the slit. The results 
of the wavelength calibration are confirmed by comparing 
the reduced Th-Ar spectra with an artificial image based 
on the known wavelengths of the Th-Ar emission lines (cf. 
figure 4). The typical calibration error of the spectra is 
less than 1 pixel, corresponding to 5 A or 1.2 A in the 
low- or high-resolution mode, respectively. 

3.2. Initial background subtraction 

Since the usual FMOS observations are carried out us- 
ing an ABAB nodding pattern of the telescope. A — B 



Fig. 4. Combined ID image of the Th-Ar spectra (at bot- 
tom) compared with an artiflcial image based on the known 
wavelengths of the Th-Ar emission lines (at top). The agree- 
ment in the position of the lines shows that the conversion 
from the column numbers to wavelengths is correct. 

simple sky subtractions can be performed using two dif- 
ferent sky images: A^ — Bn-i and An — Bn (where An 
denotes the image taken at position A of the n-th pair). 
For IRS2, these sky-subtracted frames are calculated by 
the ramp fitting algorithm applied to the sky-subtracted 
sub- frames. When the brightness of the OH airglow varies 
monotonically, most of it can be canceled by merging these 
two images according to An — B = {An — Bn)w -\- {An — 
Bn-i){l — w). The weight w is chosen to make the sum of 
the absolute count \An — 5| a minimum within the range 
— 1 < < 2. Figure 5 shows a pair of sky-subtracted im- 
ages and the merged one. Typically, the weight w is equal 
to about 0.5. 



3.3. Corrections of detector cross talk, bias difference, 
and bad pixels 

After the initial background subtraction, cross talk is 
removed by subtracting 0.15% for IRSl and 1% for IRS2 
from each quadrant. Next, the bias difference between the 
quadrants (as indicated in figure 5, the lower left quadrant 
tends to show a higher bias level) is corrected to make the 
average over each quadrant equal. After flat fielding us- 
ing the detector flat image, the registered and temporally 
prominent bad pixels are rejected, together with the four 
adjacent pixels by interpolating the surrounding pixels. 

3.4' Distortion correction and residual sky subtraction 

The processed image is converted into a combined ID 
image based on the distortion parameters obtained in the 
preparatory reduction process. However, only one line of 
each spectrum (9 pixels in width) is extracted instead of 
summing the 9 lines. One can thereby make a set of 9 
combined ID images each of which consists of a different 
part of each spectrum (cf. figure 6). In other words, the 
PSF of a fiber is divided into 9 pieces, only one of which 
is used in a combined ID image. The flexure or tempera- 
ture change in the spectrograph causes a small difference 
between the dome-flat and the scientific images along the 
vertical direction in position. This difference is corrected 
using the vertical position of the spectra of bright stars. 
In this way, the counts of the residual sky becomes smooth 
along the columns after the relative throughput correction 
of the fibers. The residual airglow lines are fitted and sub- 
tracted in these images, and then the relative throughput 
difference is multiplied to restore the noise level back to 
the original state. The residual subtracted images are re- 
combined to form an image in which the PSF of the fiber 
is determined. Medium-level bad pixels having smaller 
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Fig. 2. Correction of the spatial distortion of the dome-flat image. The y-axis of the left-hand image is converted to straighten the 
spectra, as shov^n at right. The thick lines indicate the shape of the spectra at the top and bottom edges of the detector. The left 
and right half regions are separately corrected in this lov^-resolution mode, because the corresponding mask mirrors are separated 
(Kimura et al. 2010). The projected images of each half region along the x-axis are shown along the sides of the image. 




Fig. 5. Canceling the residual OH-airglov^ emission. Left panel: Sky- subtracted image using the previous image as the sky image 
(An — Bn-i)- Center panel: Sky-subtracted image using the next image as the sky image (An —Bn)- Right panel: Merged image 
(An — B) with a weight oi w = 0.472 in which the negative spectra cannot be used because they have merged information with 
various weights. 
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Fig. 8. Relative dome-flat image to correct the throughput 
differences between spectra. Top panel: Flux-normalized 
dome-flat spectra in which the throughput differences of 
the flber have been removed. Bottom panel: Relative 
dome-flat image created by normalizing the top image along 
the columns. This image is displayed within a range of 0.8 
(black) to 1.2 (white) to emphasize the relative differences. 

size than the PSF of a fiber are replaced at the end of this 
process. 

3.5. Combine and residual background subtraction 

There are two ways to observe scientific targets with 
the ABAB nodding pattern: "normal beam switching" 
(NBS), in which all targets are observed at position A (all 
the fibers are supposed to observe "blank" sky at position 
5), and "cross beam switching" (CBS), in which less than 
half of the fibers are allocated to the targets at position A 
while the others are at position B. In the CBS observation 
mode, the same targets appear in both images collected in 
positions A or B. However, the target spectra in position 
B are merged to minimize the absolute fiux of the residual 
airglow lines in the initial background subtraction process. 
Thus the same reduction process has to be applied to the 
images in position 5, replacing the negative spectra of 
position A with the corresponding parts of the position 
B image (as shown in figure 7). The merged images (or 
all of the position A images in the NBS mode) are then 
combined into one averaged image. 

Finally, the combined image is divided into a set of 9 
combined ID images again, in order to perform a fine 
subtraction of the residual background. The recomposed 
image also goes through a bad pixel rejection process. 

3.6. Mask edge correction and CBS combine 

Although the correction for spectral distortions was per- 
formed to straighten the images of the airglow masks (as 
in figure 3), there remain small differences among the 
throughput patterns of the spectra because of the local 
shape errors of the mask elements or because of the pres- 
ence of dust on the mask mirror. These differences are 
corrected by dividing the averaged image by a "relative" 
dome-fiat image in which the common spectral features 
have been removed by normalizing the count along lines 
and columns (cf. figure 8). After correcting the relative 
differences of the throughput patterns among the spec- 
tra, the corresponding spectra from positions A (positive) 
and B (negative) taken in the CBS observation mode are 
ready to be combined. The negative spectra in the image 
are then inverted and rearranged so that they can be com- 
bined with the corresponding positive spectra (as shown 
in figure 9). 



3.7. Object mask 

In the process of fitting and subtract the residual back- 
ground, some fraction of the object fiux can be subtracted. 
The best way to retain the object fiux is to mask the ob- 
jects during the fit. Faint objects to be masked are there- 
fore selected by a human eye on the combined image and 
the reduction process is repeated with these masks ap- 
plied. The mask density should not be too high to make 
the residual background subtraction work well: the max- 
imum density allowed is roughly 75% of the image. 

3.8. Square-noise map 

Before the mask edge correction process, the distribu- 
tion of noise in an image is truncated Poisson distribution 
caused by the bad pixel rejection process in three sub- 
sequent stages. Although the noise level in an image is 
almost homogeneous at this stage, the noise level map be- 
comes complicated during the throughput correction and 
CBS combine. To estimate the noise level of each pixel, 
a frame is defined that consists of the squared noise level 
measured along columns just before the mask edge cor- 
rection. Here, the noise level is measured by a 3cr clipping 
algorithm iterated ten times for each column, while some 
extra contribution proportional to the count is added to 
the clipped pixels. This square- noise frame is divided by 
the squared image of the relative dome fiat, and reduced 
to half wherever a pair of spectra is averaged during the 
CBS combine. Note that the noise level of the bright ob- 
jects is probably underestimated because the systematic 
uncertainty (i.e. sub-pixel shifts of the spectra in raw im- 
ages and tiny variation of the throughput pattern caused 
by the instrument status such as temperature) is the ma- 
jor contribution for them. 

4. Flux calibration and check process 

4-1' Outline 

To calibrate the fiux of the scientific targets, at least one 
bright {J{AB) = 15-18mag) star in each science frames is 
needed as a spectral reference, because it must have al- 
most the same atmospheric absorption feature as that for 
the scientific targets in the same field of view. All the 
spectra are divided by the reference spectrum, and then 
multiplied by the expected spectrum of the reference star 
whose fiux and spectral type are known or can be deter- 
mined by the observed values in two different wavelength 
regions. If the cataloged or estimated spectrum of the ref- 
erence star is correct, all the observed spectra will then be 
calibrated accurately. In the next subsection, the method 
to estimate the reference star spectrum is described. 

4-2. Template stellar spectra 

Since the fiux and slope of the spectrum of a reference 
star are determined from the measured counts in an im- 
age, one needs the slope-removed template spectra includ- 
ing the intrinsic absorption features of the star. We an- 
alyzed 128 stellar spectra in the IRTF Spectral Library^ 



http://irtfweb.ifa.hawaii.edu/spex/IRTF_SpectraLLibrary/References.html 
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Fig. 6. Residual background subtraction and recomposition process. Left panel: A set of 9 combined ID images from different part 
of the PSF. Center panel: Residual subtracted images by performing fits along the columns in each image. Right panel: Recomposed 
image from the 9 residual subtracted images. 




Fig. 7. Merging the CBS images from positions A and B. Left panel: Positive spectra in the reduced An — B image in which the 
removed negative (position B) spectra contain merged information (cf. figure 5). Center panel: Negative spectra in the reduced 
A — Bn image in which the removed positive (position A) spectra are merged one. Right panel: Merged An — Bn image. 




Fig. 9. Combining a pair of target spectra taken in the CBS observation mode. Left panel: Final A — B image in which the 
positions of the targets observed only at position B are masked. Center panel: Inverted and rearranged negative spectrum that will 
be combined with the corresponding positive spectrum. Right panel: Combined spectrum. The S/N ratios of the positive spectra 
are 1.4 times higher than the ratio of the corresponding negative spectra in the combined image. 
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(Rayner et al. 2009) split into seven groups: F0-F9IV/V 
(19 objects), G0-G8IV/V (14 objects), K0-K7IV/V (10 
objects), F0-F9I/II/III (21 objects), G0-G9I/II/III (31 
objects), K0-K3I/II/III (17 objects), and K4-K7I/II/III 
(16 objects). In each group, the spectra were averaged to 
improve the S/N ratio and divided by the results of the 
linear fitting, so that slope-removed template spectra of 
seven different types were prepared (cf. figure 10). Here, 
only F, G, and K type stars were used to make the slope- 
removed template spectra, because 1) their slope-removed 
spectra become roughly straight in A - plot, 2) these 
stars are quite popular and easy to select in the target 
field, and 3) A and earlier type spectra are neither avail- 
able in the IRTF Spectral Library, nor in other database 
with similar qualities. 

Next, the correlations between the spectral types and 
the slopes of the spectra have to be established. Figure 
11 shows the distribution of measured value of 128 stellar 
spectra from figure 10. The stellar types are numbered 
from (FO) to 29 (K9) while the slopes are defined by 



slope = 



F,(1.55)-F,(1.21) 



(1) 



(1.55-1.21)F^(1.31) 

in /im~^. The correlations between the spectral types and 
the slopes are determined by second-order polynomial fits 
to the distributions of stellar types III and V: 



slopeiu = 0mi32type^ 
slopey = 0.000456t?/pe^ 



-0.0229type- 
-OmbStype- 



0.666, 
-0.654. 



(2) 



As a result, any spectrum from FO to K9 can be syn- 
thesized by multiplying a linear spectrum having a de- 
fined slope with the intrinsic absorption spectrum from 
the interpolation of the nearest two slope-removed tem- 
plate spectra. 

4.3. Calibration process 

The first step in the calibration process is the rela- 
tive throughput correction of fibers in which the averaged 
throughput of position A and B is used for the merged 
spectra in the CBS combine process. The next step is 
to remap the pixels with an increment of 5 A/pixel (1.25 
A/pixel in the high-resolution mode) based on the cor- 
relation between the observed wavelength and the pixels 
determined in the preparatory reduction process. In this 
remapping process, the observed wavelengths are multi- 
plied by the count of each pixel in order to convert the 
value from photon count to Fj^. Next, the converted 
count is divided by the atmospheric transmittance func- 
tion of the airmass present during the observation, so as 
to roughly correct for the effects of atmospheric absorp- 
tion. After this correction, the flux of the reference star is 
estimated from the converted count at around 1.31/im, as- 
suming that the total system efficiency under good seeing 
condition is 2.5% including losses at the entrance of the 
fibers. The slope of the spectrum is then measured with 
a fixed efficiency ratio between 1.21 and 1.55 /im. (The 
value of this ratio will be confirmed in the check process, 
along with the total system efficiency of 2.5%.) Finally, 
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Fig. 10. Stellar spectra in the IRTF Spectral Library. The 
slope-removed averaged spectrum of each stellar type (thick 
black lines) is used as the intrinsic absorption template. 




Fig. 11. Correlation between a stellar type and the slope of 
its spectrum. The stellar types are numbered from (FO) to 
29 (K9), and the slope is defined by equation (1). The open 
symbols are measured from figure 10: type I (triangles), type 
II (diamonds), type III (squares), type IV (pentagons), and 
type V (circles). The distributions of types III and V are fit to 
second-order polynomials represented by the thick solid and 
dashed lines, respectively. The thick dotted line indicate the 
intermediate correlation between types III and V used when 
the stellar type of a reference star is unknown. 

the observed scientific spectra are divided by the reference 
spectrum and multiphed by the stehar spectra from the 
measured flux and slope. Here, the type V and III ab- 
sorption templates are adopted as the reference spectra, 
respectively bluer than G5 and redder than Kl (cf. fig- 
ure 11). As a consequence, two of the intrinsic absorption 
templates of FV, GV, KlIII, and K4III (in figure 10) are 
used to interpolate the absorption of the reference spec- 
trum. If the stellar type of the reference star is known, 
the corresponding synthesized spectrum is used instead of 
this predicted spectrum. 

A resulting wavelength- and flux-calibrated image 
is shown in figure 12. The size of this image is 
1800x1800 pixels, the wavelength in /im is A = 0.9 + 
0.0005 X {C olumnN umber — 1) , the n-th spectrum is lo- 
cated between ^ = 9 x {LineN umber) — 8 and 7/ = 9 x 
{LineN umber) ^ and the count is in /iJy. The square- noise 
frame is converted in a similar way except that all factors 
are multiplied twice. The ID spectrum of each object is 
extracted from this image with a user-defined mask of 9 
pixels, together with the square- noise frame. An example 
of the final ID spectrum is shown in figure 13. 

4-4' Check process 

The results are checked by comparing the resulting flux 
in the J and H bands with the photometric data in the 
catalog. If all the factors contributing to the system effi- 
ciency are normal, the flux should be consistent with the 



Fig. 13. Example of a processed ID spectrum. The spec- 
trum and the corresponding la noise level are represented by 
the thick continuous line and the thin dotted line, respectively. 
The wavelength has been converted to the rest wavelength us- 
ing the cataloged redshift value (z = 0.844). 

catalog value. The accuracy of the efficiency ratio be- 
tween 1.21 and 1.55 /im can also be confirmed by these 
diagrams for most cases. Figure 14 shows an example of 
such a comparison. However, weather conditions, tele- 
scope and instrument focus, extended (non-stellar) mor- 
phology of targets, and position accuracy of the catalog 
may cause the lower observed flux of some targets than 
that expected from the catalog. Under sufficiently good 
conditions, the observed flux should match the catalog 
magnitude as shown in the thick line in figure 14 with 
some downward scatter due to small allocation error of 
the fibers. Also, if the estimated efficiency ratio between 
1.21 and 1.55 is not correct, the points of J and H will 
have small offset in this figure. 

5. Summary and Conclusions 

The reduction process for FMOS images includes two 
special processing steps. One is the segmented process- 
ing of the spectra to handle a given part of the PSF 
as a unit, while the other is the automatic modeling of 
the reference spectrum to calibrate the scientific targets. 
The segmented processing enable to keep the original 2- 
dimensional information which has a large effect on bad 
pixel filtering and detection of faint emission-lines. Most 
of the processes are carried out automatically, but the 
object mask preparation and the reference star selection 
require user judgement. The FIBRE-pac is available from 
the FMOS instrument page^ of the Subaru web site, to- 
gether with the sample dataset presented in this paper. 
The base reduction platform is IRAF and the reduction 

http:/ /www. naoj.org/Observing/Instruments/FMOS/ 



FMOS image-based reduction package 
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Fig. 12. Calibrated spectral image (at left) and a square-noise frame (at right) with a size of 1800x1800 pixels. The wavelength 
range is 0.9-1.8 /xm with an increment of 5 A/pixel. Each of the 200 spectra has a width of 9 pixels. 

scripts and sources are free and open, so that users can 
check what is happening in each step by sending the com- 
mands one at a time. 
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Fig. 14. Comparison of the observed flux with the catalog 
magnitude. The triangles and circles represent the values in 
the J and H band, respectively. The thick solid lines indicate 
the ideal values without loss of flux. 



